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The diflfusivity of Ga and Al adatoms on the (2x4) reconstructed GaAs(OOl) sur- 
face are evaluated using detailed ab initio total energy calculations of the potential 
energy surface together with transition state theory. A strong diffusion anisotropy 
is found, with the direction of fastest diffusion being parallel to the surface As- 
dimer orientation. In contrast to previous calculations we identify a short-bridge 
position between the two As atoms of a surface dimer as the adsorption site for Al 
and Ga adatoms. 



1 Introduction 

Heteroepitaxy of GaAs and AlAs on (001) oriented substrates has been suc- 
cessfully applied to create structures of reduced dimensionality. In spite of the 
considerable progress in experimental studies the knowledge of the underlying 
microscopic processes is rather poor. For example, even for a key process like 
the Ga adatom diffusimi on GaAs(OOl), the experimentally derived values for 
the activation energyldo scatter significantly (between 1.1 and 4.0 eV) and the 
very basic question "which is the direction of fastest diffusion on GaAs(OOl)?" 
has not yet been answered, unambiguously. A similar lack of understanding 
exists with respect to the question of how the diffusion properties of Al and 
Ga adatoms differ. This is crucial for controlling the growth of high quality 
GaAs/ AlAs quantum devices. 

In this paper we present calculations of the diffusivities (activation energies 
as well as pre-exponential factors) of Ga and Al adatoms on the flat GaAs(OOl) 
surface using density functional total energy calculations. 

2 Method 



Typically the rate constant for adatom hops between neighboring sites is by 
orders of magnitude smaller than the highest phonon frequencies. Therefore, 
an ab initio molecular dynamics simulation is not a practicable approach to 
evaluate the diffusivity. It is clearly more appropriate to evaluate the Born- 
Oppenheimer potential energy surface (PES) and to analyse its properties. Our 
study of the diffusivity proceeds in three steps. At first we determine the PES 
using density functional theory and the generalized gradient approximation B 
for the exchange and correlation potential. This calculation identifies |t|he local 
minima {A/} of the PES. We then apply transition state theory □ (TST) 
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to determine the individual jump rates Tj j between different sites A/ and 
A J. In a third step we use a random walk model which assumes uncorrelated 
jumps between the sites {A/} to obtain the difFusivities in the two orthogonal 
directions, [110] and [TlO]. These are the main axes of the diffusion tensor 
according to the symmetry of the considered surface structure. 
In TST the jump rate between two sites A/ and Aj is 

Tij = rljexp{-AE/kBT) , (1) 

where AE is the difference between the energies of the adatom at the site A/ 
and at the saddle point between the sites A/ and Aj. To find these sites and 
the energy barrier AE, we calculate the PES for a diffusing adatom: 

i?PES(Xi,ri)= min E{{R,}) . (2) 

^l'{^i}i=2..-.N 

The set {Ri} describes the positions of the atoms where Xi, Yi and Zi are the 
coordinates of the adatom parallel and perpendicular the surface, respectively. 
i?({Ri}) is the total energy of the surface with an adatom. 

The pre-exponential factor j is determined within the harmonic ap- 
proximation of TSTB: 



ro = 1 det\Ki\ 

V^^^^^ det\K's\ ■ 

Kj and Kg are the force constant matrices for the adatom at the site A/ 
and at the saddle point between A/ and Aj. Both matrices contain the force 
constants for Ri and for the coordinates Ri of substrate atoms next to the 
adatom. According to the TST, Kg does not contain the mode perpendicular 
to the transition plane between / and J. 



3 Results 



We assume in our study that the surface is reconstructed with the (2x4) f32 
structure which consists of rows of two As dimers and two missing dimers in 
the topmost layer. This structure is found to be present under typical growth 
conditions u and is also a structure of lowest surface energy. B 

Our calculations show that for both adatom species Al and Ga there exist 
two sheets of the PES which are relevant for diffusion. The occurrence of 
several sheets is not too surprising because in the combined configuration space 
of the adatom and the substrate there may be several saddle points between 
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two local minima A/ and Aj. Accordingly the PES as defined in Eq. (|^) does 
not need to be a unique function. 

The two sheets of the PES differ in the way the adatom interacts with 
the surface As-dimers. Either the surface displays only small relaxations (in 
particular the dimer bonds remain intact) when the adatom diffuses across the 
surface or the adatom breaks a surface dimer and causes a strong relaxation. 
For most positions {Xi,Yi) of the adatom the two PES are identical. The 
differences occur only around the the surface dimers and are caused by an 
energy barrier to break the dimers. 

If the adatom does not break the surface dimers (see Fig. |l|(a)) the long- 
bridge position Ai between surface dim£rs is the favoured site for adsorption. 
This agrees with previous calculations □ which found a long-bridge position 
as the adsorption site. However, we find that the site of lowest total energy, 
i.e., the actual adsorption site, is the short-bridge position A4 (Fig. [|(b)). 
If the adatom is adsorbed at this position the surface dimer is broken and 
the adatom forms directional bonds with the As atoms of the broken dimer. 
Such directional bonds are not formed when the adatom is located at one of 
the long-bridge sites Ai or A2 where the adatom does not break any surface 
bond. Therefore, the adatom can interact at an Ai or A2 site only with the 
completely filled dangling bonds of the As-dimers. Accordingly the adsorption 
energies at the long-bridge sites are by about 0.6 eV smaller than at the short- 
bridge sites (A3, A4). 

In order to obtain the diffusivity the jump rates on both sheets of the 
PES between the various sites of local total energy minima were combined by 
using a random walk model on a periodic lattice with inequivalent sitesu The 
resulting mobilities along the [110] and [TlO] direction are D^^^^ — 0.015 cm^/s 
exp(-1.2eV/fcBT ) and Djno] =0.03cmVs exp(-1.5eV/fcBT) for Ga adatoms 
and O.OBcmVs exp(-1.4eV/A:BT) and O.OScm^/s exp(-1.6eV//cBT) for AI, 
respectively. These diffusion constants are significantly smaller than jprevi- 
ously published values determined by ab initio total energy calculationsB This 
difference is due to the fact that the activation energy of diffusion is mainly 
determined by the energy needed to jump away from the adsorption site, which 
was not correctly identified in former investigations. 

The calculated diffusion constants exhibit a pronounced anisotropy. As it 
is already obvious in Fig. ^ the barriers for diffusion between two adjacent 
A4 sites are smaller along the [110] direction than along the [110] direction. 
Accordingly, diffusion is fastest along the [TlO] direction. Comparing the PES 
of Ga and AI adatoms we find no qualitative differences. But the PES for an 
AI adatom is more strongly corrugated and therefore it has higher diffusion 
barriers. This is a consequence of the larger adatom-As bond strength for AI 
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Figure 1: Potential energy surfaces for a Ga adatom on the GaAs(OOl)— (2x4) j32 surface, 
(a) PES sheet obtained when the As surface dimers remain intact, (b) PES sheet obtained 
when the surface dimers are broken by the adatom. The contour— line spacing is 0.3 eV. The 
atomic positions of the clean surface are indicated for atoms of the upper two layers and for 
the As-dimers in the third layer (As: filled circles, Ga: empty circles). 



as reflected by the higher cohesive energy of AlAs compared to that of GaAs. 
Therefore, AI diffuses slower than Ga in spite of the fact that AI has a higher 
pre-exponential factor (see Eq. (^)), which is mainly caused by the lower mass 
of AI as compared to Ga. 
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